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We present and analyze new exact gyraton solutions of algebraic type II on a background which 
is a static, cylindrically symmetric Melvin universe of type D. For a vanishing electromagnetic 
field it reduces to previously studied gyratons on Minkowski background. We demonstrate that 
the solutions are a member of a more general family of the Kundt spacetimes. We show that 
the Einstein equations reduce to a set of mostly linear equations on a transverse 2-space and we 
discuss the properties of polynomial scalar curvature invariants which are generally non-constant 
but unaffected by the presence of gyratons. 
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I. INTRODUCTION 

Gyraton solutions represent the gravitational field of a 
localized matter source with an intrinsic rotation which 
is moving at the speed of light. Such an idealized ul- 
trarelativistic source can describe a pulse of a spinning 
radiation beam and it is accompanied by a sandwich or 
impulsive gravitational wave. 

The gravitational fields generated by (nonrotating) 
light pulses and beams were already studied by Tolman 
[ij in 1934, who obtained the corresponding solution in 
the linear approximation of the Einstein theory. Ex- 
act solutions of the Einstein-Maxwell equations for such 
'pencils of light' were found and analyzed by Peres 
and Bonnor [J-Q. These solutions belong to a general 
family of pp -waves d, 0] ■ 

In the impulsive limit (i.e., for an infinitely thin beam, 
and for the delta-type distribution of the light-pulse in 
time) , the simplest of these solutions represents the well- 
known Aichelburg-Sexl metric Q which describes the 
field of a pointlike null particle. Subsequently, more gen- 
eral impulsive waves were found [ol-fisj (for recent reviews 
see [ll[l3). 

The gyraton solutions are generalization of -waves 
which belong to the Kundt class for which the source — 
the beam of radiation — carries not only energy, but also 
an additional angular momentum. Such spacetimes were 
first considered by Bonnor in [l8| , who studied the grav- 
itational field created by a spinning null fluid. In some 
cases, this may be interpreted as a massless neutrino 
field [lil. 

Gyratons on Minkowski background are locally isomet- 
ric to standard -waves in the exterior vacuum region, 
outside the source. The interior region contains a nonex- 
panding null matter which possesses an intrinsic spin. In 
general, these solutions are obtained b y k eeping nondiag- 
onal terms gui in the Brinkmann form (20| of the pp -wave 
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solution, where u is the null coordinate and are or- 
thogonal spatial coordinates. The corresponding energy- 
momentum tensor thus also contains an extra nondiago- 
nal term T!„j — ji. In four dimensions, the terms g^i can 
be set to zero locally, using a suitable gauge transfor- 
mation. However, they can not be globally removed be- 
cause the gauge invariant contour integral § gui{u, x^) da:* 
around the position of the gyraton is proportional to the 
nonzero angular momentum density ji, which is nonvan- 
ishing. 

These gyratons were investigated (in the linear approx- 
imation) in f2]| in higher dimensional fiat space, the ex- 
act gyraton solutions propagating in an asymptotically 
flat D-dimensional spacetime were further investigated in 
[23 | . They proved that the Einstein's equations for gyra- 
tons reduce to a set of linear equations in the Euclidean 
{D — 2)-dimensional space and showed that the gyraton 
metrics belong to a class of so called VSI spacetimes for 
which all polynomial scalar invariants, constructed from 
the curvature and its covariant derivatives, vanish identi- 
cally 23] . (For the discussion of spacetimes with nonvan- 
ishing but nonpolynomial scalar invariants of curvature, 
see [2JI-) Subsequently, charged gyratons in Minkowski 
space in any dimension were presented in [25| . 

In [2^, the exact gyraton solutions in the asymptoti- 
cally anti-de Sitter spacetime were found. Namely, they 
obtained Siklos gyratons which generalize the Siklos fam- 
ily of nonexpanding waves [27i | (investigated further in 
[2^1) which belong to the class of spacetimes with con- 
stant scalar invariants {CSI spacetimes) [29l - [3^ . 

Recently, the large class of gyratons on the direct- 
product spacetimes was found in [11], where we showed 
that this class of gyratons has similar properties as the 
previous gyratonic solutions: the Einstein's equations re- 
duce to a set of linear equations in transversal 2-space and 
these spacetimes belong to the CSI class of spacetimes. 

Let us also mention that string gyratons in supergrav- 
ity were recently found in 3j|. Supersymmetric gyraton 
solutions were also obtained in minimal gauged theory 
in five dimensions in [35i] , where the configuration repre- 
sents a generalization of the Siklos waves with a nonzero 
angular momentum in anti-de Sitter space. 
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The gyratons are important in studies of production 
of mini black holes or in cosmic ray experiments. The 
theory of high energy particle collisions was developed in 
[sgI [33] and was applied to gyraton models in [sl] . 

The main purpose of this paper is to further extend the 
family of gyratonic solutions. In particular, we present 
new gyraton solutions of algebraic type II, propagating 
in the Melvin universe. 

The Melvin universe [s^, H^l is a nonsingular electro- 
vacuum solution with physical properties which are in- 
teresting both from a classical and a quantum point of 
view. The spacetime represents a parallel bundle of mag- 
netic (or electric) flux held together by its own gravita- 
tional attraction. The transverse space orthogonal to the 
direction of the flux has a nontrivial spatial geometry. 
It was represented in [41| by a suitable embedding dia- 
gram which resembles a tall narrow-necked vase. Also 
m it was shown that no motion can get too far 

from the axis of symmetry. This aspect is analogous to 
the attractive effect of a negative cosmological constant 
in the anti-de Sitter universe. 

The Melvin universe was considered as an important 
model in astrophysical processes related to gravitational 
collapse because of its stability. It was shown in [10, Elj 
that the spacetime is surprisingly stable against small 
radial perturbations and also against large perturbations 
which are concentrated in a finite region about the axis of 
symmetry. The asymmetries are radiated away in gravi- 
tational and electromagnetic waves [i^, HJ] . 

The Melvin universe also appears as a limit in more 
complicated solutions, in [2^ it is obtained as a specific 
limit of a charged C-metric. The Melvin universe has 
been generalized to Kaluza-Klein and dilaton theories 
[4^, to nonlinear electrodynamics [i^l, and has impor- 
tant applications in the study of quantum black hole pair 
creation in a background electromagnetic field [48-54] . 

The Melvin universe recently attracted a new inter- 
est because it is possible to find gravitational waves in 
the Melvin universe [i^ by an ultrarelativistic boost of 
the Schwarzschild-Melvin black hole metric [i^j. It was 
shown that these wave solutions are straightforward im- 
pulsive limits of a more general class of Kundt spacetimes 
of type II with an arbitrary profile function, which can 
be interpreted as gravitational waves propagating on the 
Melvin spacetime. The gyraton spacetimes investigated 
in this paper are generalizations of such Kundt waves 
when their ultrarelativistic source is made of a 'spinning 
matter'. 

The paper is organized as follows. In Section |TT] we re- 
view basic information about the Melvin universe which 
will be useful in the paper. We derive the ansatz for 
the gyraton metric by a direct transformation from the 
Kundt form of the metric to Melvin's coordinates. We 
also review the transverse space geometry of the wave 
front. 

In Section IIIIl we derive the field equations and we 
simplify them introducing potentials. We discuss the 
structure of the equations and the gauge freedom of the 



solutions. 

Next, in Section IIVI we solve the Einstein-Maxwell 
equations in the special case of the 0-independent space- 
times, especially with a thin matter source localized on 
the axis of symmetry. 

In Section |V] we concentrate on the interpretation of 
the gyraton solutions. We discuss the properties of the 
scalar polynomial invariants and the geometric properties 
of the principal null congruence. We evaluate the curva- 
ture tensor in an appropriate tetrad, discuss the Petrov 
type, the matter content of the spacetime, and properties 
of the electromagnetic field. 

The main results of the paper are summarized in con- 
cluding Section IVH Some technical results needed to de- 
rive the field equations, spin coefficients and invariants 
are left to Appendices [B] and O 

II. THE GYRATONS ON THE MELVIN 
SPACETIME 

A. The Melvin universe 

In this section we briefly review basic properties of 
the Melvin spacetime [H, |4^ H^l which will be useful 
throughout the paper. The Melvin universe describes an 
axial electromagnetic field concentrating under influence 
of its self-gravity. The strength of the electromagnetic 
field is determined by the parameters E and B. In cylin- 
drical coordinates (t, z, p, (/>), the metric and the Maxwell 
tensor read 

ds^ = Y.'^i-dt^ + dz^ + dp^) + Y.-^p^d(t)'^ , (2.1) 
F = Edz Adt + BY.-'^pdphdcl) , (2.2) 

where 

S = 1 + J^?EM/3' . (2.3) 
The constant p^M is given by the parameters B as^ 
g,,, = ^{E^ + B^). (2.4) 
Introducing double null coordinates, 

we obtain an alternative expression for metric (12.11) and 
the Maxwell tensor (1^ : 

ds^ = s2(-2dudv + dp2) + S^V^d0^ , (2-6) 
F ^ Edv Adu + BY.-^pdpAd(j) . (2.7) 



a = SttG and £0 are gravitational and electromagnetic constants. 
There are two standard choices of geometrical units: the Gaus- 
sian with K = 8iT and £0 = l/47r, and SI like with x = Eo = 1- 
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The electromagnetic field can be rewritten also in the 
complex self-dual form,^ 

T = B{dvAdu- i^-^p dp A dcj)) . (2.8) 

with the complex constant B defined as 



B = E + iB 



(2.9) 



The metric (|2.1I) jesembles a vacuum solution of the 
Levi-Civita family [6] for a large value of p. For p^M = 
{E, B — 0) spacetime reduces to the Minkowski space- 
time in cylindrical coordinates. For E 0, B — the 
Maxwell tensor describes an electric field pointing along 
the z-direction, whereas ioi E = 0, B we get a purely 
magnetic field oriented along the z-direction. 

The metric admits the four Killing vectors 



dt , dz 



, zdt + tdz vdv - udu 



(2.10) 



which correspond to staticity, cylindrical symmetry, 
and invariance under a boost transformation. Us- 
ing the adapted null tetrad k = d^, I = Y^~^du, and 
m — ^{'S^^dp — i'Ep~^d^), the only non-vanishing 
components of Weyl and Ricci tensors are 



^'2 = ^Qem (-1 + ^QemP^)^' 

;r f V^— 4 



(2.11) 



This demonstrates that the Melvin universe is a non- 
vacuum solution of the Petrov type D, except at points 
satisfying p = 2/^Pem, where the Weyl tensor vanishes. 
It is interesting to note that the scalar curvature vanishes, 

To conclude, the Melvin spacetime belongs to the fam- 
ily of non-expanding, non- twisting type D electrovacuum 
solutions investigated by Plebahski [55|. As a conse- 
quence, it also belongs to the general Kundt class 
which will be important in the following text. 



B. The ansatz for the gyratons on Melvin universe 

Gyratons are generalized gravitational waves corre- 
sponding to null sources with intrinsic rotation. In gen- 
eral, the gyraton solutions are obtained by adding non- 
diagonal terms to the metric of the standard gravita- 
tional wave solutions, or in other words, by keeping the 



^ We follow the notation of 6] , namely, T = F + i *F is a com- 
plex self-dual Maxwell tensor, where the 4-dimensional Hodge 



dual is *_F, 



^i' per 



The self-dual condition reads 



-kT = —iT. The orientation of the 4-dimensional Levi-Civita 
tensor is fixed by the sign of the component £vup4> = pS^. The 
energy-momentum tensor of the electromagnetic field is given by 
T — £o F P^ 

-t — 2 M ^P' 



non-diagonal terms Qui in the standard Kundt metric [g'] . 
Therefore, we derive the ansatz for the gyraton on Melvin 
spacetime by adding such new terms to the Kundt form 
of the Melvin metric. It can be explicitly obtained from 
(|2.6[) by transformation 

v = Y.-'^r, (2.12) 



which leads to 

ds^ = -2du dr + ds^ 4- 2rWi dudx' . 
Here we introduced a 2-dimensional metric 



dsi 



^-^p^dct)^ 



and r-independent 1-form W 



Widx^ 



W 



2^dp. 



(2.13) 



(2.14) 



(2.15) 



These tensors can be understood as tensors on space 
spanned by two coordinates p, 4>. This space can be cov- 
ered by other suitable spatial coordinates x% and we will 
use the Latin indices i, j, ... to label the corresponding 
tensor components. 

By an appropriate transformation of coordinates 0, 
i44j , the 2-dimensional metric ds\ can be transformed into 
a conformally flat form, which in the standard complex 
null coordinates C, Q reads 

2 



dsi 



p2 



dCdC 



(2.16) 



Such a transformation brings the metric (j2.13l) into the 
Kundt form. 

The gyraton generalization of (I2.13|) then reads 

-2Y.'^Hdu'^ - 2dudrH 
- 2{rW,-J:'^a,) dudx' 



ds^ = 



dsi 



(2.17) 



We have added the term —21^^11 du^ which represents a 
gravitational wave on the Melvin universe H^l with 
an arbitrary profile function iJ, and the non-diagonal 
terms — 2E^aidMdx* characteristic for gyratons. It will be 
shown in the following that these terms can be generated 
by specific gyratonic matter. 

Transforming back to the Melvin coordinate v and 
cylindrical coordinates p, 4>, we obtain the ansatz for the 
metric describing the gyraton on Melvin spacetime,'^ 

ds^ = -2^^Hdu^ - 2Y?dudv 4- {Y?dp^ + S^^^d-^^N 



21]^ 



,dudp^ 



(2.18) 



The function H{u,v,p,(l)) can depend on all coor- 
dinates, but we assume that the functions ai{u,p,4>) 
are v independent (it actually follows from the Maxwell 
equations as will be shown below). Let us note that 
the previously cited works assumed also the function H 
V independent. 



^ Here we use notation different from |43l . |44|. we use 
instead of —Hdu^ to match our notation in 1331. 



-2Hdu'^ 
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C. The ansatz for the matter 

The metric should satisfy the Einstein equations with 
a stress-energy tensor generated by the electromagnetic 
field and the gyratonic source, 

G^. = >^(T™ + T-;) . (2.19) 



We assume that the electromagnetic field (|2.7|) modified 
by the gyraton is given by 

F ^Edv Adu + BY.-'^pAp hd(f) + cTjdu A dx^ . (2.20) 

Similar to (33[ , we have added the term aj du A dx^ . 

To evaluate Maxwell equations, it is useful to write 
down the self-dual form of the Maxwell tensor F. The 
Hodge dual of (|Z!20)) reads 



= S dw A d-it - EYj-^p dp A ^ 



{*a + Ba — E *a)j du A da;-' 



(2.21) 



where the star * means the 2-dimensional Hodge dual de- 
fined on the transversal space, see the next section (jll Pp . 
For the self-dual Maxwell tensor F we thus obtain 



(2.22) 



F {dv A du - iS^^p dp A dcj) 

+ {S — i *a).j du A dx^ ^ . 

where we introduced a complex transverse 1-form Sj , 

BSj = {(T + i*a)j + iB{a + i*a)j . (2.23) 

This form is self-dual with respect to the Hodge duality 
* on the transversal space, 

*Sj^-iSj, (2.24) 

and therefore it can be written using a real 1-form sj: 

Sj = {s + i *s)j . (2.25) 

The original 1-form ai can be expressed in terms of Sj as 

aj ^ Esj - B*{s- a)j . (2.26) 

In the following we use Sj(v,u, p, 4>) as a basic variable 
for the electromagnetic field. 

The stress-energy tensor T™ corresponding to the field 
([2:221) is given in (|A3|) . 

Finally, we must define the gyratonic matter by spec- 
ifying the structure of its stress-energy tensor. It is ob- 
tained from the standard stress-energy tensor of a null 
fluid by adding terms corresponding to 'internal spatial 
rotation' of the fluid: 

>f T^^'- = j„ du^ + 2jp du dp + 2j0 du dcf) . (2.27) 

We admit a general coordinate dependence of the source 
functions j„(f , u, p, if) and jj(v, u, p, 4>). However, it will 



be shown below that the field equations enforce a trivial 
V dependence. 

The gyraton source is described only on a phenomeno- 
logical level, by its stress-energy tensor (|2.27p . which is 
assumed to be given, and our aim is to determine its in- 
fluence on the metric and the electromagnetic field. How- 
ever, we have to consider that the gyraton stress-energy 
tensor is locally conserved. It means that the functions 
j„ and ji must satisfy the constraint given by 







(2.28) 



Of course, if we had considered a specific internal struc- 
ture of the gyratonic matter, the local energy-momentum 
conservation would have been a consequence of field equa- 
tions for the gyraton. Without that, we have to require 
([2:28)) explicitly. 

To conclude, the fields are characterized by functions 
S, H, aj, and Sj, which must be determined by the field 
equations provided the gyraton sources ju and jj and the 
constants E and B are given. 



D. The geometry of the transverse space 



The geometry (|2.18p identifies the null geodesic con- 
gruence generated by which is parametrized by 
an affine parameter u, the family of null hypersur- 
faces u — constant, and 2-dimensional transverse spaces 
u,v = constant. 

The gravitational wave moves along the null direction 
dy, i.e., it propagates with the speed of light along the 
z-direction, which is the direction of the electromagnetic 
field. The hypersurface u = constant corresponds to the 
surface of the constant 'phase', and the transverse spaces 
u,v — constant are spatial wave fronts of the wave. 

Physical quantities do not depend on the affine param- 
eter V, or this dependence is trivial and it will be explic- 
itly found. Specifically, the geometry of the transverse 
space is v independent. 

It turns out to be convenient to restrict various quan- 
tities to the transverse space. For example, we can inter- 
pret flj and Si as components of u-dependent 1-forms on 
the transverse space. Our goal is to formulate all equa- 
tions for physical quantities on the transverse spaces. For 
that we need to review some properties of the transverse 
geometry. It was studied in detail in [13], and on a gen- 
eral level in 33| , nevertheless it will be useful to mention 
some of the properties explicitly. 

The transverse metric is obtained by restriction of the 
full metric (|2.18p to the transverse space and it is given 
by the expression (|2.14p . 



dsl = grijdx'dx^ = S^dp^ + £-2^2^02 



(2.29) 



The associated Gauss curvature is given by the scalar 
curvature i?T, 



K 



1 



E4 



(2.30) 
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It is obvious that only the electromagnetic field is respon- 
sible for the non-flatness of the transversal space — it is 
insensitive to the presence of the gyraton. For q^m — 
the curvature vanishes and we get the flat plane. 

In general, the curvature is not constant and it is fi- 
nite everywhere. The Gauss curvature K has maximum 
on the axis p = where it is equal to 2gEM', it is pos- 
itive for < p < 2^/2/^1 



1, and vanishes on the cir- 
cle at p — 2\/2/^Pem- For p > 2\/2/ y/g^, it goes to 
negative values and it has its minimum K — at 
p = 2A/3/Y^fe^. Then it grows again, and as p — > +oo 
the curvature vanishes, K ^ Q~ . 

The circumference of a circle of constant radius p is 
vanishing when p — ^ +oo. Therefore, we measure a much 
shorter circumference for larger p — as if we would move 
"along the stem of the wine-glass toward the narrowing 

end," go,!!!,!!!. 

The 2-dimensional Levi-Civita tensor associated 
with the metric (|2.29p is e = pdpA dep. The covariant 
derivative will be denoted by a colon, e.g., at.j. We raise 
and lower the Latin indices using grij, which differs from 
lowering indices using gap thanks to non- vanishing terms 
gui- We use a shorthand for a transverse square of the 
norm of a 1-form a, as 



(2.31) 



In two dimensions, the Hodge duals of 0, 1 and 2-forms 
ip, ai, and /y , respectively, read 



1 



1 



^ ^^ij , (*a)» = aj-e^j , *f = T:fi3^'^ = -fp<l> ■ 

I p 

(2.32) 

For convenience, we also introduce an explicit notation 
for 2-dimensional divergence and rotation of a transverse 
1-form a,-, 



div a = ai 



1. ^ 



1 



ap 0,p 



rot a 



(2.33) 



For 2-form fij we get 

1 1 

div / = .Uj'-^ = ;^(/0p.p--/0p)d0H — jfp^.^dp , (2.34) 
1j p p 

and rot f — 0- We can generalize the action of divergence 
and rotation also on a scalar function / as div / = 
and rot/ = — *d/. Note that the divergence and ro- 
tation are related as div a = rot *a, and the relation to 
the transverse exterior derivative is da = *rot a. Clearly, 
div div a — 0, div rot a = 0, and rot da — 0. 
The Laplace operator of a function reads 

1 1 

(2.35) 

and for a transverse 1-form rj it is defined as 
A ?7 = d div rj — rot rot rj. 



Finally, the transverse space is topologically trivial 
since it has topology of a plane. We can thus assume that 
the Poincare lemma {duj = w = da) holds, which in 
terms of rotation and divergence means that divcj = 
implies lj — rot cr. However, since the transverse space is 
non-compact and we do not know a priori boundary con- 
ditions for various quantities at infinity p ^ 00, we have 
to admit non-trivial harmonics. Therefore, we cannot 
assume a uniqueness of the Hodge decomposition. More- 
over, in some cases it can be physically relevant to con- 
sider also topologically nontrivial harmonics which are 
singular, e.g., at the origin p = 0. Such solutions would 
correspond to fields around singular sources localized on 
the axis. However, we will ignore these cases in a general 
discussion. 



III. THE FIELD EQUATIONS 
A. The field equations for matter 

Now, we will investigate the equations for matter, i.e., 
the Maxwell equations for electromagnetic field and the 
condition ()2.28|) for the gyraton source. 

Both Maxwell equations for real Maxwell tensor are 
equivalent to the cyclic Maxwell equation for the self- 
dual Maxwell tensor ([2^ . 



dT^B 



dy (s + i *(s — a)) du A du A dx^ 
- (rot s + i div (s — a)) du A e =0 



(3.1) 



From the real part we immediately get that the 1-form Si 
is V independent, d^Si = 0, and rotation- free. 



rot s — 



(3.2) 



From the imaginary part it follows that the 1-form is 
also V independent (as we have already mentioned above) 
and it satisfies 



div (s — a) = 



(3.3) 



Equations p.2p and (13. 3|) guarantee the existence and 
determine the structure of potentials which will be dis- 
cussed in detail in Section fill CI 

Next, we analyze the condition (|2.28l) for the gyraton 
source. When translated to the transverse space, it gives 

- dyj, dx' + {-dyju + div + Y^a'd^j.,) du = . 

(3.4) 

The source functions ji must be thus v independent and 
ju has to have the structure 



j„ = t;div (E^j) -I- L 



(3.5) 



where l{u,x'') is a w independent function. The gyra- 
ton source (I2.27P is therefore fully determined by three 
D-independent functions l(u,x^) and ji(u, a;^). 
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Equation (|3.5|) gives us also an insight into interpre- 
tation of the gyratonic terms jV They are composed 
from two contributions: one is related to a kind of 'heat 
flow' which changes energy ju of the fluid, and the other 
which is related to intrinsic rotation of the fluid. The 
source representing 'heat flow' has thus non- vanishing di- 
vergence div (S^j) and we require a vanishing rotational 
part rot(E^j). In opposite, the source representing in- 
trinsic rotation has vanishing heat flow, i.e., it satisfies 



div (S^ = 



(3.6) 



Such a source can be written in terms of a rotational 
potential v as j ^ — S^^rotj/. In components it means 



1 



Jp = 



S4 



(3.7) 



Physically more relevant is the rotational part of the 
source, since it can describe the spin of the null fluid, 
or, in a specific limit, of the polarized beam of light. 
Terms related to heating flow have bad causal behav- 
ior and therefore typically do not satisfy various energy 
conditions and they are thus rather unphysical. Inter- 
pretation of the gyraton source was discussed previously 
also in [2l|,[2l,[26[|3l■ 



B. The Einstein equations 



The Einstein gravitational law (|2.19l) needs the Ein- 
stein tensor and the electromagnetic and gyraton stress- 
energy tensors. These quantities can be found in the 
Appendix]^ We can combine them and inspect various 
components of the Einstein equations. 

The wM-component determines the function S, namely 
it gives the condition 



P(S,p)^ 



2si; p = 



(3.8) 



It is straightforward to check that it is satisfied again by 
S in the Melvin form (lOI) . 

The transverse diagonal components pp and (fxp require 



where we introduce the exterior derivative fij of the 
transverse 1-form a; and its Hodge dual b{u,x^), 

fij = - = (* b)^j , (3.12) 
b = *f^ rota . (3.13) 

In terms of b, equation (|3.11|) can be rewritten as 

j = ^ rot (S^ b) + E^dg + 2g^^{s - a) . (3.14) 

Here and in the following dg represents just the trans- 
verse gradient dg = g^i dx* , and for simplicity we skipped 
the transverse indices. 

It is useful to split the equation into divergence and 
rotation parts by applying div and rot : 



div(E2 j) = div(s2dg). 



(3.15) 



rot (S^j) = --A (E^'fe) + rot (E^dg) - 2feM b . (3.16) 



Here we have used the relations p.2|) and (|3.13|) . The 
formula (I3.15P is the equation for g, p.l6p is the equation 
for b and together with (I3.13P it determines a^. In the 
next section we will return to these equations introducing 
suitable potentials which allow us to escape the necessity 
of taking an additional derivative of (I3.14[) when deriving 
the equation for a^. 

Finally, from the uu-component of the Einstein equa- 
tion we obtain 

ju = div i^^dg) v + Y?{^h- {j:-^),php) 

+ -J^'^b^ + 2Y?a'g., + a„div (E^a) -I- 5 div (E^a) 



SgEM (s - af 



(3.17) 



When we compare the coefficient in front of v with p.lSp 
we find that it has structure consistent with (13. 5p . The 
nontrivial ^-independent part of (|3.17l) gives the equation 
for the metric function h, 

V V ),p ,P! 2 (3.18) 

— 9udiv (E^a) — g div (E'^a) + 2gEM (s — o)^ . 



diH = 



(3.9) 



C. Introducing potentials 



thus we obtain the explicit v dependence of the metric 
function H as 

H = gv + h, (3.10) 

where we have introduced ^-independent functions 
g{u,x^) and h{u,x^). 

The remaining nontrivial components of the Einstein 
equations are those involving the gyraton source (|2.27l) . 
The uz-components give the equation related to ji, 

h = U'f^^'-{^'),kg'''fJ^+g., + ^{s^-a,) , (3.11) 



In the previous section we have found that the Maxwell 
and Einstein equations reduce to two potential equations 
p.2|) . p.3|) . and two source equations p. lip . p.l7|) . 

According to the two dimensional Hodge decomposi- 
tion we can express the 1-form using two scalar poten- 
tials k{u,x^) and X{u,x^), 

a = dK-|-rotA. (3.19) 

These potentials control the divergence and the rotation 
of tti as 

diva = AK, rota = -A A. (3.20) 
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Comparing with (I3.13|) we thus obtain the equation for 
A in terms of 6, 

AA = -6. (3.21) 

The first potential equation (|3.2I) gives immediately 
that Si has a potential (p{u,x^), 

s^dif . (3.22) 

Equation (j3.3|) implies that there exists a potential 
'ip{u,x^) satisfying 

s — a = — rotV' ■ (3.23) 

In terms of these potentials the 1-form (j2.26|) from the 
real Maxwell tensor (|2.7p reads 

cr = £;d</j + BdtP . (3.24) 

The potential and i/' are not, however, independent. 
Substituting and into ((X^ we obtain the 

key relation among the potentials: 

d{ip - k) + rot (t/- - A) = . (3.25) 

If the Hodge decomposition was unique, the gradient and 
rotational parts would be vanishing separately, i.e., we 
would get ip = K and = X (up to unimportant con- 
stants). The non-uniqueness of the Hodge decomposi- 
tion is linked to the possible existence of a non-trivial 
harmonic 1-form r]i{u,x^), 

Ari = 0, (3.26) 

in terms of which the gradient and rotation parts of 
(|3.25p can be expressed as 

if — K = divTy , (3.27) 
V'-A = -rot?7. (3.28) 

These are equations for electromagnetic potentials (p 
(or -0, respectively) in terms of the metric potentials k 
and A. The 1-form rj encodes an extra freedom, which 
allows a nontrivial electromagnetic field not uniquely de- 
termined by the metric. (Such contributions would allow 
one to take into account, for example, an additional elec- 
tromagnetic charge localized at the origin of the trans- 
verse space, cf. the discussion of particular cases in sec- 
tion |lVl) However, sufficiently restrictive conditions at 
the infinity and the smoothness on the whole transverse 
space for the potentials would eliminate this freedom, so 
the case 77 = is a rather representative choice. 

After eliminating the electromagnetic potentials, we 
need to formulate the equations for k and A. We start 
with the divergence part of p.l4p which can be rewritten 
using the modified Laplace operator acting on the metric 
function g: 

S]2(Ag-(S-2)^^5^^) ^div(E2j) . (3.29) 



However, for g solving this equation, (I3.15|) is also the 
integrability condition for the quantity T,^{dg — j). It 
can thus be written in terms of a potential w, 

S^(dg- j) =rotw , (3.30) 

or, more explicitly, 

dw = -E2(rot.g + *j) . (3.31) 

For the source p.7p without intrinsic 'heating', the 
right-hand side of p.29p is zero and the function uj has 
an additive contribution from the rotational potential v, 
namely it has to satisfy 

d(a; - i^) = -E^rotg . (3.32) 

The function lu contains information from the source 
J and from the metric function g relevant for the rota- 
tional part of equation (|3.14l) . Indeed, substituting the 
potentials and rot uj into (j3.14l) we obtain 

rot (^^"6 - 2gEMV' + = . (3.33) 

Substituting (|3.21|) . p.28p . and integrating (absorbing an 
integration constant to lo), we derive the equation for the 
potential A 

is'' A A + 2giEMA = w + 2gEM rot 77 . (3.34) 

Taking into account relations (I3.28|) and (|3.26p , we obtain 
the alternative equation for -0 

iEVi^ + 2&M^ = cj. (3.35) 

The other metric potential n remains unrestricted. 
This non-uniqueness is related to the gauge freedom dis- 
cussed in detail in section UlI El This coordinate freedom 
allows us to set the potential n to an arbitrary convenient 
form, e.g., to eliminate it completely. 

D. Discussion of the field equations 

We have thus formulated all field equations as equa- 
tions on the transverse space. They are written in a sep- 
arated form, i.e., they can be solved one after the other: 
First, one has to find harmonic 1-form rj satisfying equa- 
tion p.26p and metric function g satisfying p.l5p . It 
allows one to integrate the function lo which together 
with rot rj appears as a source in equation p.34p for the 
potential A. Using the gauge freedom one can choose the 
other potential k. Equations p. 271) and p.28p determine 
electromagnetic potentials and through equation (I2.26P 
the electromagnetic field. Finally, the remaining metric 
function h is determined by equation p. 181) , in which the 
previously computed quantities contribute to the source 
on the right-hand side. 
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The most complicated field equations — (|3.29|) . p. 181) . 
and p.34p — are partial differential equations on the 
2-dimensional space, which are solvable, at least in princi- 
ple. They all contain a modified Laplace operator, equa- 
tions (|3.29p and p.lSp for g and h the same one, namely 

I]2(a / - /^^) = div (S^d/) . (3.36) 

Solutions for particular cases (assuming rotational sym- 
metry) will be discussed in section IIVI 

It is important to observe that, except equation (13.18^ 
for h, the field equations are linear. We can thus super- 
pose two solutions simply by adding the fields together. 
Only in the last step, when computing the source for 
equation p.lSp . one has to include total superposition of 
the fields g, a^, and Ui since the expression for the source 
is non-linear. 

Finally, we have not paid much attention to the 
u dependence of the studied quantities. All metric func- 
tions, matter fields and sources can depend on the co- 
ordinate u and this dependence does not enter the field 
equations except in one term on the right-hand side of 
equation p. 181) . The profile of the gyraton in the u di- 
rection can thus be specified arbitrarily. It corresponds 
to the fact that both matter and gravitational field move 
with the speed of light and information on different hy- 
persurfaces u — constant evolves rather independently. 

Also the dependence of the metric and fields on the 
coordinate v is very simple and it was found for all quan- 
tities explicitly. 

E. The gauge transformation 

The coordinate transformation v ^ v — v — xi^iX-') 
accompanied by the following redefinition of the metric 
functions and matter fields: 

9 = 9, h = h + gx + duX , at = hi - x,% , 

■^i X.i I "'i ^ X-i J ,„ 

(3.37) 

K = fi-X, A = A, ip^ip-x, V' = V' , 
uJ = oj , V = V , 
jt=ji, t^t + xdivj, 

leaves the metric, the Maxwell tensor, and the gyraton 
stress-energy tensor in the same form. Therefore, all the 
field equations remain the same. This transformation is 
thus a pure gauge transformation and we can use it to 
simplify the solution of the equations. 

There are two natural choices of gauge: we can elimi- 
nate either the metric potential k or the electromagnetic 
potential f. In the first case a — rot A and (p — div rj. 
In the latter case s = 0, a = Bd4>, k — —divr/, and 
a = rot ip. 

Let us mention that in [3^ an analogous gauge trans- 
formation allowed us to choose also the metric function 



g. For the gyraton on the Melvin universe the metric 
function g decouples from n and it is gauge independent. 

The discussed gauge transformation has a clear geo- 
metrical meaning: it corresponds to a shift of the origin 
of the affine parameter v of the null congruence dy. Note 
that such a change redefines transverse spaces. 



IV. SPECIAL CASES 

In this section we will study the special solutions of 
the field equations. Namely, we restrict to the axially 
symmetric situation, i.e., to the case when the geometry 
and the fields are invariant under action generated by 
the rotational vector d^. Further, we concentrate on the 
gyraton generated by a thin beam of matter concentrated 
at the origin of the transverse space which means on the 
axis of symmetry. 

Thank to linearity mentioned at the end of section 
IIIIDi we can discuss various special cases separately. 
However, the geometry of the spacetime with gyraton 
is not merely a superposition of the individual contribu- 
tion since the nonlinear coupling in the metric function 
h. 

We do not discuss in detail the u dependence of the 
fields. It does not enter the field equations, except in the 
source term of equation p.l8p for h. The u dependence 
of the gyraton sources and corresponding dependence of 
other fields can thus be chosen arbitrarily. 

The symmetry assumption enforces that quantities a^, 
g, h, Ui, ji, and l are ^-independent. It can induce a 
slightly weaker condition on the potentials: typically, 
they have a linear dependence on 0. 

The thin beam approximation requires that gyratonic 
matter is concentrated at the origin o of the transverse 
space given by p = 0, i.e., ji and l should be distribu- 
tions with the support at the origin. However, we relax 
this condition slightly in the case of gyratonic 'heat fiow' 
discussed in section HVBl 

In all discussed cases we use a natural gauge 

K = , i.e., a = rot A . (4-1) 

Together with the symmetry assumptions it implies 
A,0 = constant. 



A. Pure gravitational gyraton 

We start with the simplest vacuum case: we set ji — 
and t = and we assume no pure electromagnetic contri- 
bution, i.e., 77 = 0. The equation for g has only a trivial 
solution 9 = 9o = constant, the function uj is also a triv- 
ial constant and we obtain equation p.34p with vanishing 
right-hand side. Taking into account that A^^ = constant 
we obtain that A must be (f) independent (a (/)-linear term 
in A would require an analogous term in the source) and 
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we obtain an ordinary differential equation 

4(?E 



l( P 
P 



(4.2) 



Substituting (|2.3p for S, it is possible to obtain two in- 
dependent solutions, one regular at the origin, 



A = -7I]-l(l-|^^™p^) , 



(4.3) 



and the other behaving as log p near the origin, which 
corresponds to a delta source at the origin and it will be 
discussed in section ITVDI 

The metric l-form a = rot A has thus components 



and its 'strength' b = rot a is then 



= -7^5- ( 1 - 40™P 



(4.4) 



(4.5) 



These can be plugged into (|3.18l) which turns out to be 

^(p^p),p = -7^^(i - y..p') (1 - i&mp^ 
^ . . . ^^-^^ 

It can be integrated explicitly; however, the result is 
rather long, so we skip it. 



B. Non-spinning beam with 'heat' flow 



As we discussed in section lTlI Al the gyraton source can 
have two contributions: one, which changes gyr atonic en- 
ergy density ju and the other corresponding to intrinsic 
rotation. Let us investigate the case when the gyratonic 
energy is concentrated at the origin (a thin beam) but 
there is axially symmetric energy flow in the dp direc- 
tion which accumulates energy at the beam (a 'heating' 
process). Namely, we assume that 9„j„ = div (E^j) is 
nonzero only at the origin, so elsewhere the transverse 
flow ji must satisfy equation p.6p . Such ji has the form 



a I . 



(4.7) 



and the increasing energy of the gyraton is given by 



ju ^ av So 



(4.8) 



Here, So stands for the transverse delta function localized 
at the origin o, normalized to the standard metric volume 
element on the transverse space. 

Since the 'heating' is localized only at the origin, the 
gyraton source (|4.7p can be locally written using the 
source potential i^, cf. (13.71) . as 



a 



(4.9) 



Note however, that the potential cannot be defined glob- 
ally and it is not well behaved at the origin. 

We again assume no pure electromagnetic contribu- 
tion, i.e., rj = 0. The requirement of the axial symmetry 
enforces that the difference E~^rotw between dg and j, 
cf. p.30p . must be zero, i.e., cj = 0. The 0- independent 
metric function g must thus satisfy g,p = jp which gives 



ZTT 



(4.10) 



The equation for A has again the form (|4.2p and in 
this case we choose the trivial solution A = 0. With the 
gauge K = we thus obtain the only nontrivial field to be 
the metric function g, otherwise = and ai = 0. The 
source for equation p.lSp is also trivial, given only by i, 
and wc will study it in the next case. 



C. Non-spinning light beam 

A particular example of the gyraton source is standard 
null fluid. The thin non-spinning beam localized at the 
origin is described by the source 







(4.11) 



We can set all the fields except h to be zero: g = 0, 
ai = 0, and <7i — 0. The equation for h outside the origin 
is p~^{ph,p),p = which (with proper fixing of the source 
constant) gives 



h^—\ogp. (4.12) 



D. Thin gyraton — spinning light beam 

Finally we proceed with the most characteristic repre- 
sentant of the gyratonic matter. It is a simple null beam 
of energy localized at the origin with no heating, which, 
however, contains an intrinsic energy rotation. Since we 
have a point-like source at the transverse space, we can 
speak about inner spin instead of a global rotational en- 
ergy flow. 

The gyraton source has a form'^ 



j ~ —I] ^rot ly 



L^edo, (4.13) 
, with V = . (4.14) 



The symmetry assumptions together with the no-heating 
requirement implies g — constant, which we choose to be 
zero in this case. Equation p.30p also implies that uj = v. 



^ The exact structure of the singular source ji at the origin can be 
read out from the singular solution l|4.16|l of equation l|3.34ll for 
A below. 
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Ignoring again the electromagnetic contribution, = 0, 
we obtain the equation for A 

The solution of the homogeneous equation with a singu- 
lar behavior ~ log p corresponding to the delta function 
(I4.14[) at the origin reads 

+ i(l-|foMP^)log(ifoMP^)) . (4-16) 

The multiplicative constant in the argument of the loga- 
rithm can be chosen arbitrary since it generates only an 
additional homogeneous contribution of the form (j4.3p . 
The metric 1-form and its rotation b are 

Op = , 

l3 ^ / 324 136 148 

- feMP^lOgd&MP^)) , (4.17) 

and 

b = 2^feMS-5((l-|eEM/5^) log(if?EMP^) 

2324 13 61 
+ 2 + QemP — -^QemP ~ ssI^'emP ) ■ 

The source for equation p.lSp becomes cumbersome and 
lengthy, but treatable, in principle. 



(/50, ipcf, being constants and (p and functions of p only. 
Substituting into p.25p we get 

(f.P ~ y ^0) dp + (</50 + ^i^,p) d0 = , (4.20) 
which implies 

V,P^—'4''P, V'.p^-— (4-21) 

Taking into account (|2.3p . it can be easily integrated and 
substituting back to (|4.19p we obtain 

tp = X - ip^{\ogp+ IqemP'^) +il^^(t) ■ 

At this moment it is easier to solve equation p.35p 
for -0 with vanishing right-hand side instead of equation 
([334)1 for A. It has solutions in the form (|4?3| and (|4?T6| . 
The metric potential A is then given by the second of the 
equations (|4.19p . The metric potential k is vanishing 
thanks to our gauge. 

After choosing the solution for tp, we can thus com- 
pute all quantities a, b, s, and a and substitute them to 
equation (|3.18l) for h. 

Let us mention that solution (|4.19p is singular at the 
origin. A careful distributional calculation would show 
that equations for the potentials (I3.22p and (|3.23p may 
not be satisfied at the origin. Tracing this singular term 
back to Maxwell equations, it could lead to non-vanishing 
electric charges localized at the origin. However, since 
we have not written down the Maxwell equations with 
sources, we do not discuss these terms in more detail. 



E. Electromagnetic wave 

In the previous examples we have ignored the possibil- 
ity of a nontrivial electromagnetic field. Namely, we have 
assumed that the electromagnetic field is given by the 
metric potentials via relations tp = k and ^p = \. How- 
ever, we have already observed that equation (|3.25p ad- 
mits also other solutions which we have parametrized us- 
ing a harmonic 1-form 77. To include such solutions we 
should classify all 1-form harmonics on the transverse 
space. But if we restrict to the axially symmetric fields 
we can solve the potential equation p.25p directly, with- 
out referring to 77 explicitly. 

Let us study a pure electromagnetic contribution to 
the matter, i.e., we assume t = 0, ji = here. We can 
thus take a trivial vanishing solution for g which implies 
w = 0. 

The symmetry assumptions tell us that 1-forms s, a, 
and a are (j) independent, which implies that all potentials 
including ip — k and ijj — X can be at most linear in </): 

(p- K = ifip) + (f^(j> , t/j - X^ tp{p) +ij^(j> , (4.19) 



V. PROPERTIES OF THE GYRATON 
SPACETIMES 

A. Gravitational field 

In this section we discuss some of the geometrical prop- 
erties of the gyr atonic solutions. 

One of the important characteristics of spacetimes are 
the scalar polynomial invariants which are constructed 
only from the curvature and its covariant derivatives. It 
was shown that gyratons in the Minkowski spacetime 
(2^ have all the scalar polynomial invariants vanishing 
(VSI spacetimes) [23^1, the gyratons in the anti-de Sitter 
[26| and direct product spacetimes (33)] have all invariants 
constant (CSI spacetimes) [29| . 

In these cases, the invariants are independent of all 
metric functions which characterize the gyraton, and 
have the same values as the corresponding invariants of 
the background spacetime. We observe that similar prop- 
erty is valid also for the gyraton on Melvin spacetime, 
however, in this case the invariants are generally non- 
constant, namely, they depend on the coordinate p. This 
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property is a consequence of the general theorem hold- 
ing for the relevant subclass of the Kundt solution, see 
Theorem II. 7 in [sgI ). 

Values of some of the scalar curvature invariants can 
be found in Appendix [Cl 

The metric ()2.18p admits the null vector k ^ dv It is 
a Killing vector for g — and div (S^j) = 0, i.e., for the 
no 'heating' part in the gyratonic source. The covariant 
derivative of k is given by 

k^,p = -^-^{d,H)k^kfi + E-2fc[„V^]E2 . (5.1) 

We observe that the congruence is not even recurrent 
as in the case of a gyraton on direct product spacetimes 
[ssj . For = we recover the formula in Garfinkle 
and Melvin }43| . In general, the non-recurrency of the 
congruence is related to the non- vanishing spin coefficient 
Tnp = — j^S^p, cf. d^l, calculated in Appendix IB] 
The null character of k and the condition (15.11) im- 



ply that the null congruence with tangent vector k is 
geodesic, expansion-free, sheer-free and twist-free, and 
the spacetime thus belongs to the Kundt class. 

Next we calculate components of the curvature ten- 
sors with respect to the following adapted null tetrad 
{k, I, m, m} p] 



k — , 



m 



(5.2) 



^/2S 



do 



P 



Here, we have introduced the projection of a transverse 
1-form a on the vector m 

am = ra^ai = ap — i — — (a + i *a)p , (5.3) 

and we will use an analogous notation also for compo- 
nents of the transverse gradient of a real function / 



/,m — ^ I ,i f ,p i f,<p 1 /,m — f.r. 



(5.4) 



The dual tetrad of 1-forms {e^^), Q'-^\ Q'-"'\ e*")} has 
a simple form 

e^'^' ^dv + Hdu - a , e(') = SMu , 



Hdu- 
- i *dp] 



(5.5) 



Q(m) ^ (dp + i*dp) 

v2 



With respect to this tetrad, we have found that the 
non-vanishing curvature components are given by four 
new components and by those which are the same for the 
Melvin universe (|2.6p . The non- vanishing Ricci scalars 
are 

1 1 



*12 = 



4^/2 S3 
1 1 
2S2 



^22 = ^^{AH-{j:-')^pHp 



2,2 



(5.6) 



2a*5,, + S-2(^ + 5^)div(I]2a) 



and the non- vanishing Weyl scalars read 
1 



aV2 S3 



1 1 

'^A — 

^ 2S5 
+ S(gH 



2S Qr, 



du) 



■g.r. 



'■p,p 



S( H^pp~—^H^rt> 

I — (a„ 
P 



S2 

-2i — H , 
P 



S4 



H,p + duttp + gap 



(5.7) 



-I- 2i — {H^^+dua^+ga^] 

In particular, there exists a relation between ^'3 and $12: 
1 1 



$1 



^-3 = 



4V2 S3 



[-i(S2),p5 + 4.g,wi] 



(5.8) 



Therefore, the metric (|2.18p describes the transversal 
gravitational wave (^"4 term) in the k direction with a 
longitudal wave component (vl/3 term) . The gravitational 
wave is accompanied by an aligned pure radiation field 
($22 term) with non-null component ($12 term) propa- 
gating in the Melvin universe. In fact, the scalars ^'3 and 
<I'i2 are generated by the gyratonic functions ai and the 
function g. In general the spacetime (j2.18p is of Petrov 
type //. 

Let us now investigate the subcases of our solutions. 
When we set the gyratonic functions = 0, the Ricci 
scalars become 

1 1 



$12 



$22 = 



2V2S3 
1 1 



5,1 



2EV2 

and the Weyl scalars then read 
1 1 



[pipHp)^p + E^H,^^] 



(5.9) 



<I»3 = 



2V2S3 

'^4 = l^A^pH,pp 



2E5p 
1 1 
2^ 



+ (2pS,p - S)iJ,p 



(5.10) 



1 1 



2S3p^ 



We again obtain the spacetime with similar characteris- 
tics as for the full gyratonic metric (|2.18p . However, the 
scalars 5*3 and $12, which now depend only on function 
g, are now related by an even simpler relation: 

¥12 = *3. (5.11) 

If we assume additionally H to be w-independent (i.e., 
g = 0) we obtain the only non- vanishing scalars $22 and 
^'4 in the same form as in (|5.9I) and (|5.10p . This case 
and its subcases were thoroughly discussed in ^44i].^ 



The terms $22 and ^4 have a httle different form which is caused 
by the slightly different choice of null tetrad in our paper. 
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Finally, let us mention that for Pem = the background 
reduces to Minkowski spacetime and we thus recover the 
gyraton moving on the Minkowski background as an im- 
portant subcase of our solutions. 

B. Electromagnetic field 

The gyraton propagates in a non-null electromagnetic 
field (|2.22p . the influence of which on the geometry is 
characterized by its density ^em 

((^ . The electro- 
magnetic field is modified by the gyraton through the 
CTidu A dx' terms. The electromagnetic field can be 
rewritten in terms of potentials using p.24p as 

F = E(dv Adu + duA dw) 

(5 12) 

+ B{^-^pdpAd(j) + duAdij). 

It describes a superposition of electric and magnetic 
fields, both pointing along the z direction, which are 
modified by the gravitation field of the gyraton. The ad- 
ditional term does not have a simple structure of electric 
or magnetic field, however both are of the form du A df 
with / being the proper potential. 

The electromagnetic field projected on the null tetrad 
(|5.2p is characterized by three scalars $i, 

(5.13) 

It follows that the non-null electromagnetic field is 
aligned with the principal null direction k of the grav- 
itation field, but this vector is not a double degenerate 
vector of the field. 



VI. CONCLUSION 

We have derived and analyzed new gyraton solutions 
moving with the speed of light on electro- vacuum Melvin 
background spacetime in four dimensions. This solution 
extends the gyraton solutions previously known on the 
Nariai, anti-Nariai, and Plebahski-Hacyan universes of 
type D, and on conformally fiat Bertotti-Robinson and 
Minkowski space. 

The gyraton solutions describe a gravitational field 
created by a stress-energy tensor of a spinning (circu- 
larly polarized) high-frequency beam of electromagnetic 
radiation, neutrino, or any other massless fields. The 
gyratons generalize standard gravitational -waves or 
Kundt waves by admitting a non-zero angular momen- 
tum of the source. The interpretation is that the null 
matter in the interior of the gyratonic source possesses 
an intrinsic spin (or non-zero angular momentum). This 
leads to other nontrivial components of the Einstein 
equations, namely, Gui in addition to the pure radia- 
tion MM-component which appears for pp -waves or Kundt 
waves. 



We have shown that it is possible to define the gyra- 
ton by adding the gyratonic terms Ui to the gravitational 
wave on the Melvin spacetime in cylindrical coordinates 
in a similar way as we have defined them in the gen- 
eral Kundt class. We were able to find an ansatz for the 
gyraton metric on the Melvin spacetime by direct trans- 
formation from the Kundt class of metrics (jll B[) . 

We have further demonstrated that the Einstein- 
Maxwell equations reduce to the set of linear equations on 
the 2-dimensional transverse spacetime which has a non- 
trivial geometry given by the transverse metric. These 
equations can be solved exactly for any distribution of 
the matter sources. In general, the problem has been 
thus reduced to a construction of scalar Green functions 
for certain differential operators on the transverse space. 

We have solved and analyzed the field equations for 
particular examples with the axial symmetry. In these 
cases the equations reduce to ordinary differential equa- 
tions. 

We have analyzed geometric properties of the principal 
null congruence and we have found that it is not recurrent 
contrary to the case of gyratons on direct product space- 
times. We have explicitly calculated the curvature tensor 
and determined that the gyratons on Melvin spacetime 
are of Petrov type II and belong to the Kundt family of 
shear-free and twist-free nonexpanding spacetimes. The 
gyratonic term generates the non-trivial Ricci $12 and 
Weyl V^s scalars, in addition to the gravitational waves 
investigated in [4^ . We found also a very simple relation 
(15. 8|) between these components. By studying particular 
subclasses we have shown that our solutions are general- 
izations of those from [13] . 

The scalar polynomial invariants of the metric (I2.18P 
are in general non-constant (although, some of them are 
zero) — they depend on the coordinate p. The invariants 
are not affected by the presence of the gyratons, they are 
the same as for the Melvin background. The same prop- 
erty was proved for gyratons on backgrounds belonging 
to VSI or CSI families of spacetimes. 

It would be interesting to investigate a generalization 
of our ansatz for more complicated spacetimes which 
could allow, e.g., an inclusion of a cosmological constant. 



Acknowledgments 

We wish to thank to Tomas Pechacek and Otakar 
Svitek for helpful discussions and to Marcello Ortaggio 
for his paper about gravitational waves in the Melvin uni- 
verse which motivated our work. H. K. was supported 
by Grants No. GACR-202/09/H033, No. GAUK 12209, 
and Project No. SVV 261301 of the Charles Univer- 
sity in Prague. P. K. was supported by Grant No. 
GACR 202/09/0772, and both authors thank the Project 
No. LC06014 of the Center of Theoretical Astrophysics. 



13 



Appendix A: The Einstein equations 

Here we present quantities needed for evaluation of the 
Einstein equations. 

The inverse to the metric (I2.18P is 



1 2 

^ P ^ (Al) 

1 Y? H 

+ 2{ap—dp + a^—d^)dv + {2— + a^)dvdv. 

The stress-energy tensor T"^^' of the electromagnetic 
field (|2.22p can be defined as 



(A2) 



where F = F + i-kF is the complex self-dual Maxwell 
tensor. The 4-dimensional Hodge dual is defined by 
*Fpi, = ^EpypaF'"^ and the Maxwell tensor Fp,i, satis- 
fies the self-duality condition -kj- = —iJ-. 

The non-vanishing components of the stress-energy 
tensor (IA2I) are 



S2 



H 



= 2feM + - «) ) : 
^T-;^ = ^{ap-2sp), 

pp ~ ~ ~^9pp J 



(A3) 



Here we have used only the metric (|2.18l) . without any 
usage of the field equations. 



Appendix B: The NP formalism 



Calculating the Newman-Penrose spin coefficients 
with respect to the tetrad (|5.2p , we recover again that the 
congruence k is nonexpanding and nontwisting (p^p — 0), 
sheer-free ((T^p = 0), geodesic and affine parameterized 
(knp = Enp = 0). In addition, the tetrad is gauge invari- 
ant and it is not parallelly transported along the null con- 
gruence because it does not satisfy n^p — tt^f = Enp — 0- 

The remaining spin coefficients are 



Anp — , /^Np — 2 ^ ' 
1 1 



4S2 

^NP ^ ^^{{g + du)arn + gm} , (Bl) 
_ 1 1 _ 1 1 

11, , 11 



2V2S2p^ 



2V2Sp 



rpEM 



S6 S4 



where the density ^em was defined in (j2.4p . 

The non-vanishing components of the stress-energy 
tensor (|2.27p of the gyratonic matter are 



>(T^u = Ju ^ V div (E^j) + L , 

X^ ui — • 

The Einstein tensor for the metric (|2.18p reads 



(A4) 



1 



Guu - §2^62 + S2(A H + ^^^H.p) + Y:'{dlH)a 



2/a2 r^^„2 



+ 2Y?a'dyH^i -f {d^H + du) div {T?a) 
+ 2HGuv J 

G„p = i— 5,^ - ap(G„„ - d^H) + d^H^p , (A5) 
2 p 

1 41] 

Gpp Guv 7 

= +d^H) . 



Appendix C: Scalar polynomial curvature invariants 

As we have already mentioned, the scalar curvature 
invariants are independent of all metric functions which 
characterize the gyraton, and have the same values as the 
corresponding invariants of the Melvin universe (cf. [sg]). 
Let us stress here, however, that the invariants are gen- 
erally non- constant, namely, they depend on the coordi- 
nate p. In this appendix we list some of the curvature 
invariants. 

The scalar curvature for the whole gyraton metric 
(|2.18p is zero, i? = 0. Next, we define the following 
scalar polynomial invariants constructed from the Rie- 
mann tensor: 



n(2) Dab T}cd 

n(3) _ nab ncd ne/ 

— cd^^ ef^^ ab ; 

/?(4) _ nab ncd ne/ npq 

n(5) _ nab ncd ne/ rypq nrs 
IX — il, ca-'i' ef^^ pqi^ rs^^ ab 



(CI) 



Using the GRtensor package in Maple, we get the explicit 
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2g^ 



expressions: 



EM/3 2 4 
8 e EM P 



E8 



4£L 



:m / 1 
2 Vl£ 

\^256 



^^EMP'' 



16 CBM 



3 2 4 

4 ^emP 



T^emP^ - 20 



.4 8 _ J. 3 6 
kemP ifit^EMr 



8 t^EMP 



SSpem/O + 65 



(C2) 



EM f 3 5 10 

t^EMF 



256 t^EMA 



64 cemP 



-f + 64-127fe 



8 t^EMP 



- 204 



We explicitly observe that these invariants do not de- 
pend on any of the metric functions ai, g, and h which 
characterize the gyraton. 

The invariants (|C2p mimic the behavior of the Gauss 
curvature of the transverse space discussed in detail in 
(|IID[) . They have their maximum on the axis p = 
and they are vanishing as "the neck of the vase closes off 



asymptotically" as p tends to infinity. For Pem = we get 
the identically vanishing invariants, i.e., the invariants for 
the gyratons on Minkowski background (VSI). 

In Maple tensor package GRtensor there is defined a 
set of curvature invariants CMinvars. For completeness, 
we present the explicit expressions for them: 

R^R2^ Wll ^ W2I = , 



i?3 ^ 4^ , W2R 



4E1 



28^1 



{QemP^ - 4)^ , 
j{QE»p^~Af , 



Mil = M2I = MA = MM = 



MIR= ^™ "2 



22E1 



M2R = A/3 = 



24S16 



(fiiEMP - 4) , 
(PemP^ - 4)2 , 



26E2 



(£»EM/5^ - 4)^ 



(C3) 
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